A Monte Carlo simulation study on the wetting behavior of water on graphite surface 
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This paper is an expanded edition of the rapid communication published several years ago by the 
author (Phys. Rev. B, v76, 0J^H02(R), 2007) on the simulation of wetting transition of water on 
graphite, aiming to provide more details on the methodology, parameters, and results of the study 
which might be of interest to certain readers. We calculate adsorption isotherms of water on graphite 
using grand canonical Monte Carlo simulations combined with multiple histogram reweighting, based 
on the empirical potentials of SPC/E for water, the 10-4-3 van der Waals model, and a recently 
developed induction and multipolar potential for water and graphite. Our results show that wetting 
transition of water on graphite occurs at 475-480 K, and the prewetting critical temperature lies in 
the range of 505-510 K. The calculated wetting transition temperature agrees quantitatively with 
a previously predicted value using a simple model. The observation of the coexistence of stable 
and metastable states at temperatures between the wetting transition temperature and prewetting 
critical temperature indicates that the transition is first order. 

PACS numbers: 68.35.Rh, 64.70.Fx, 82.20.Wt 
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I. INTRODUCTION 



When a fluid adsorbs on a solid surface at temperatures 
below its liquid- vapor critical temperature (T c ), the ad- 
sorbed film either spreads across the surface (wetting) or 
beads up as a droplet (nonwetting) as the pressure ap- 
proaches the saturated vapor pressure P svp of the fluid. 
Wetting transition describes the transition between those 
two kinds of behavior. Physically the wetting transition 
corresponds to the phenomena when the contact angle 
of the liquid drop on the surface changes from a nonzero 
value to zero. Analysis of wetting transition was first pre- 
sented 30 years ago by Cahn [l| and Ebner and Saam. 
Q They showed that if a fluid does not wet a particular 
surface at low temperature, then the system ought to ex- 
hibit wetting transition at some temperature T w below 
T c . In terms of adsorption isotherms, the wetting phe- 
nomenon should manifest itself as following three differ- 
ent patterns. (1) At temperatures below T w , adsorption 
beginning with a thin film increases slightly as the pres- 
sure increases towards the saturation pressure P SV p- At 
-P svp the bulk vapor condenses completely, and the ad- 
sorption coverage becomes infinite. On a coverage versus 
pressure diagram, the adsorption isotherm reaches P svp 
with a discontinuous jump (infinite slope). (2) In the 
temperature range between T w and the prewetting criti- 
cal temperature T pwc , the thin film grows as the pressure 
increases until it jumps to a thick, liquid like film of finite 
thickness at some pressure less than P svp . This thin-to- 
finitc film transition or wetting transition is followed by 
continuous growth until condensation occurs at P svp - (3) 
At temperatures higher than T pwc , the film grows con- 
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FIG. 1: Schematic diagram of adsorption isotherms near a 
wetting transition. T pwc > T 3 > T 2 > T w > Ti. 



tinuously and the prewetting transition disappears. A 
schematic diagram showing these three types of adsorp- 
tion patterns is given in FIG. [U 

Since the first theory on wetting transition was de- 
veloped a variety of experimental [3T-[lq] and theoretical 
studies fl7i - [29| have been performed on the wetting tran- 
sition of fluids on various solid surfaces. Most of these 
studies were focused on simple fluids such as He and H 2 
isotopes on alkali metal surfaces. One common feature 
of these systems is that the fluid-surface interaction is 
only weakly attractive. This implies that the bead-up 
of fluid on the surface at P svp will be favorable over the 
continuous growth of a film. 

Finn and Monson [l8| were among the first to cal- 
culate the wetting temperature of fluid on solid sur- 
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face using molecular simulations. They predicted the 
wetting behavior of Ar on solid CO2 surface using 
isobaric-isothermal Monte Carlo simulations. Shi et 
al. (30j reevaluated this system using grand canoni- 
cal Monte Carlo (GCMC) simulations plus multiple his- 
togram reweighting techniques. Errington [H( studied 
the same system by employing a new simulation method. 
Bojan et al. [25[ studied wetting behavior of Ne on sur- 
faces with various interaction strengths. Curtarolo et al. 
[26j used GCMC simulations to study the wetting be- 
havior of inert gases on alkali and Mg surfaces. Shi et al. 
[28j studied the wetting transition of hydrogen isotopes 
on Rb surface by including the quantum effects of the flu- 
ids using path integral hybrid Monte Carlo simulations. 

The wetting transition of fluid-fluid systems were stud- 
ied by both experiments and theory, [32J but the wetting 
behavior of fluids on solid surfaces was only reported 
for atomistic molecules such as inert gases. No wetting 
transition has ever been seen for any molecular fluids on 
solids other than hydrogen and its isotopes, which are 
essentially spherical molecules. To our knowledge, the 
wetting transition involving water has not been stud- 
ied until very recently 29] although water is an exten- 
sively studied molecule. It is known that water does not 
wet many surfaces (such as graphite) at room temper- 
ature. Theoretically, the wetting transition of water on 
graphite is expected to occur at a temperature below its 
bulk critical temperature. In a recent paper by Gatica 
et al., [29l | wetting temperatures of water on graphite 
had been predicted using a simple model based on the 
water-solid interaction calculated from empirical poten- 
tials. The recommended T w from their calculations for 
water on graphite is 474 K. 

In this study, we report evidence for the first-order 
wetting transition of water on graphite from molecular 
simulations. We estimate the wetting transition temper- 
ature and prewetting critical temperature of water on 
graphite using grand canonical Monte Carlo simulations. 
The paper is organized as the following: The next section 
describes the potential models and simulation method- 
ology. Section III presents the results and discussion. 
Section IV summarizes our findings. 



II. POTENTIALS AND METHODS 

Water- water interaction is described by the SPC/E 
model. [HI This model is widely used in modeling sys- 
tems involving water. The model includes a Lennard- 
Jones site located on the oxygen atom and three partial 
charge sites on each atom. The parameters for this model 
are given in TABLE fl] The critical temperature of water 
calculated using the SPC/E model is 635 K. [H| This 
value is the closest to the experimental value (647 K) 
compared with the predictions by many other popular 
nonpolarizable water potentials. [35j Up to date there 
is no potential available being able to reproduce in ev- 
ery detail the properties of real water. [35| It is known 



that the ability of an interaction potential to describe 
the bulk critical behavior is a necessary (though not suf- 
ficient) requirement in order for it to predict the wetting 
transition behavior of the fluid on a surface. [32[ There- 
fore, we chose the SPC/E model among tens of water 
potentials available. 

The graphite surface is modeled as a smooth basal 
plane. The Lennard- Jones interaction between a water 
molecule and the graphite surface is given by the 10-4-3 
potential (3|| 



V s lLj(z) = 27T£ sf of f Ap s 



V z 



10 



z 



3A(0.61A + z) s 



(1) 



where z is the distance between the oxygen atom in a wa- 
ter molecule and the graphite surface in the surface nor- 
mal direction, a s is the number density of carbon atoms 
in graphite, and A is the distance between the graphene 
sheets in graphite. The graphite surface corrugation is 
not included in this potential. The graphite-water in- 
teraction parameters £ s f and <r s f are calculated from the 
Lorentz-Bcrthelot rules, 

e sf = (e s e f ) 1/2 , er sf = (cr s + cr f )/2. 

The values of the parameters are: p s =114 nm -3 , 
A=0.335 nm, £ s =0. 05569 kcal/mol, and er s =0.340 nm for 
graphite, £f = £o=0.1554 kcal/mol, and Cf = cro=0.3165 
nm for water. 

A recently developed effective potential for the dipole- 
induced dipole, dipole-quadrupole, and quadrupole- 
quadrupole interactions between polar fluids and 
graphite (3tJ is used to calculate the water-graphite polar 
interactions. For water/graphite system, only the induc- 
tion and dipole-quadrupole terms are important, 
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where £0 is the vacuum permittivity, fee is the Boltz- 
mann's constant, fj,{ is the dipole moment of the water 
molecule, ac is the isotropic polarizability of a carbon 
atom in graphite, and 0c is the permanent quadrupole 
moment on each carbon atom in graphite. The values for 
these parameters are /Xf=1.85 Dcbye, ac = 1-76 x 10 -3 
nm 3 , [H| and 6 C = -3.03 x 10- 40 C m 2 . [H We note 
that /if calculated for SPC/E model is 2.07 Debye, but 
1.85 Debye is chosen here so that comparison can be made 
between this study and Gatica et a/.'s prediction. 

Ewald summations were applied in simulations to ac- 
count for the long-range correction to electrostatic inter- 
actions. Since only two dimensional periodic boundary 
conditions along x and y directions were applied in the 
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TABLE I: Potential parameters for SPC/E model. [33] r H is 
the O-H bond length, #hoh is the H-O-H bond angle, eo and 
o"o are Lennard- Jones parameters for the O atom, qo and qu 
are the partial charges on O and H atoms, respectively. 



roH H 


#HOH 


eo[kcal/mol] 


(TO H 


qo[e] qn[e\ 


0.1 


109.47° 


0.1554 


0.3165 


-0.8476 0.4238 



adsorption simulations, a pscudo-two-dimcnsional Ewald 
summation method [4(| was used. The total electrostatic 
energy is given by 
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where Vq is the volume of the simulation cell and a is 
the Ewald convergence parameter. In Eq. (|3]), the first 
term is the reciprocal sum for all the Gaussian charges, 
k is the reciprocal lattice vector, qj is the partial charge 
on interacting site j, r n j is the vector from r„ to Vj. 
The second term is the standard real-space sum, erfc is 
the complementary error function. The third term is 
the self-exclusion on each molecule, plus the 1-2 and 1-3 
intramolecular exclusions, where site is the total number 
of partial charge sites in the molecule, and list contains all 
the 1-2, 1-3 exclusions. The last term is the 2D correction 
term, where M , is the z component of the total dipole 
moment of the simulation box. 

We have used grand canonical Monte Carlo simulations 
combined with multiple histogram reweighting (MHR) 
method [U, |42j to compute the saturated coexistence 
chemical potentials [43| and adsorption isotherms [3(| of 
water on graphite at various temperatures. More details 
on the MHR method can be found from the original lit- 
erature fill 421 and several other articles on its applica- 
tions. [43T44| The GCMC cell for adsorption simulations 
is a rectangular box with volume of 2000(7f . The height 
of adsorption box in the z direction, H, is 15<Tf , the sides 
in x and y directions are 11.5<7f. The lower plane normal 
to the z axis is modeled as the graphite surface. 

Special care was taken to avoid capillary condensa- 
tion effects on the wetting transition properties of the 
system. The plane opposite to the graphite surface is 
modeled as a hard repulsive wall. The hard wall is al- 
ways dry to the liquid phase and wet to the gas phase. 
Therefore it helps to suppress the capillary condensa- 
tion. 0, l45T - |4a ] The capillary condensation can also 
be eliminated by using a simulation cell with sufficient 
separations between the adsorbent surface and the op- 
posite hard wall. [IH, [2^, [26|, [28| However, an overall 



small cell volume is preferred for the MHR technique. 
Therefore, we have performed series of trial simulations 
using various cell heights to search for an appropriate 
value of H. Trial simulations were performed using cell 
heights ranging from lOcxf to 40(7f at interested temper- 
atures and chemical potentials (pressures). It was found 
that the adsorption properties such as isotherms obtained 
at 15<7f are consistent with those obtained at 20<7f, 30crf , 
and 40<7f within the statistical fluctuations (see FIG. [2] 
for an example). This indicates that 15c7f is adequate for 
effectively avoiding the influence of capillary condensa- 
tion for the systems of interest. Based on this we chose 
15<7f as the cell height in most of our GCMC simulations. 
Additional simulations with if =20<7f were performed as 
verifications at each temperature. 

The type of move to attempt during a GCMC simula- 
tion was selected randomly with probability of 0.45, 0.45, 
0.05, and 0.05 for displacements, rotations, creations, and 
deletions of a water molecule, respectively. Each simula- 
tion included equilibration of 80 x 10 6 MC moves and pro- 
duction of 20 x 10 6 MC moves. Histograms were collected 
every 20 MC moves during the production. The cut- 
off of pair- wise Lennard- Jones interaction between water 
molecules was 0.9 nm as suggested by the ori gin al litera- 
ture, without long-range correction applied. [33| 

In the light of Gatica et a/.'s prediction, [29| we chose to 
perform GCMC simulations at 460, 470, 480, 490, 500, 
510 K with varying reduced chemical potentials to ob- 
tain the histograms for bulk water and water/graphite 
adsorption systems. The values of bulk saturation chem- 
ical potentials (/ig V p) a ^ each different temperature can 
be determined using the MHR method by combining the 
histograms collected. The MHR provide very precise val- 
ues of /i svp through the equal area criterion, [43| which 
is very important for studying wetting transitions. Suffi- 
cient overlap between histograms of adjacent state points 
is necessary in order to use the MHR technique. We em- 
ployed the method proposed by Shi et al. [3(| to check 
the overlap of any two adjacent state points. According 
to this method, the grand canonical partition functions 
extrapolated by histogram reweighting method for any 
two adjacent state points should approximately satisfy 
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where the subscript HR indicates that the partition func- 
tion in the numerator has been extrapolated from the 
histogram reweighting of the state point in the denomi- 
nator. The recommended value for S is 0.65 for checking 
the overlap of two adjacent state points. j3Cj Whenever 
the overlap criterion defined in Eq. ((4]) is not satisfied by 
any of two adjacent state points, additional simulations 
at state points that bridge them were performed. From 
preliminary simulations we found that the wetting tran- 
sition of water on graphite occurs approximately between 
470 and 480 K and the wetting critical temperature is be- 
tween 500 and 510 K. In order to narrow down the values 
of T w and T pwc , we performed additional simulations at 
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FIG. 2: Typical adsorption isotherms of water on graphite 
from GCMC simulations and MHR, where p* is the reduced 
number density and p* is the reduced chemical potential. The 
curves correspond to T—475, 490, 510 K, from left to right, 
which were computed from MHR using histograms collected 
from simulations with cell height H—15at. The symbols are 
data from individual GCMC simulations and were not in- 
cluded in the MHR calculations. Circles: H=15a{, squares: 
_ff=20<7f, diamonds: H=30at, triangles: H—40a{. 



475 K and 505 K. 



III. RESULTS AND DISCUSSIONS 

GCMC adsorption simulations were carried out for re- 
duced chemical potentials up to saturation under each 
selected temperature. Isotherms were calculated using 
MHR based on the histograms collected during the sim- 
ulations. Three representative isotherms are plotted in 
FIG. [2] Some of the data points calculated directly from 
GCMC simulations but not included in the MHR cal- 
culation are also plotted in the figure to compare with 
isotherms from MHR. It can be seen that the difference 
between the densities obtained from direction GCMC 
simulations and those from MHR are small. This serves 
as a test of the accuracy of the MHR isotherms. 

To present the isotherms obtained at different temper- 
atures in a concise and clear fashion, we define a param- 
eter as 



X = exp 



where fi* is the reduced chemical potential, /x* vp is the 
saturation chemical potential at the reduced temperature 
T*. Plotting the adsorption coverage versus x* gives a 
diagram similar to FIG. [TJ which helps one to identify 
the wetting transition points without ambiguity. The 
parameter x* ' 1S the ratio of the activity to the activity 
at saturation, with x* = 1 corresponding to p = /i* vp , or 
P = P svp - Additionally, x* = P/Psvp if ideal behavior is 



FIG. 3: Adsorption isotherms of water on graphite from 
GCMC simulations and MHR. The curves correspond to 
T=510, 505, 500, 490, 480, 475 K, from left to right. 



assumed in the bulk vapor phase. However, water vapor 
cannot be treated as an ideal gas under the interested 
simulation temperatures in this study. For example, the 
experimental compressib ility factor of the saturated wa- 
ter vapor is about 0.865 Hi at 500 K. 

As the chemical potential was increased toward sat- 
uration, three different types of behavior in the growth 
of the water adsorption film on graphite were observed, 
corresponding to three ranges of temperature. Adsorp- 
tion isotherms for water/graphite at several representa- 
tive temperatures are shown in FIG. [3J 

At temperatures below 475 K, the adsorption cover- 
age is minuscule until the saturation chemical poten- 
tial is reached, which indicates partial wetting or non- 
wetting. The isotherm jumps to the saturated liquid 
density at Xsvp; which can be seen from FIGs. [3J [3J 
and from the density profile growth patterns shown in 
FIG. |U The sharp increase of density between x*=0.991 
and x* — 1-008 corresponds to a first-order transition 
from nonwetting to liquid condensation (FIG. [4|. At 
X*=1.008, much of the density profile becomes compa- 
rable with the saturated liquid density profile (the sat- 
urated liquid density at 475 K is p* «0.89 from bulk 
GCMC simulations) except for the first peak corresponds 
to the density of liquid film in contact with the wall. 

In contrast, the simulation results at temperatures be- 
tween 480 K and 505 K manifest quite different behav- 
ior. Taking the isotherm at 490 K as an example, there is 
sudden jump in adsorption from minimum to a finite cov- 
erage of about p*=0.25 at x*=0.965 to x*=0.972 (FIG. 
[5]). But apparently the increased coverage does not cor- 
respond to a liquid condensation (the saturated liquid 
density at 490 K is p* w0.85). As the chemical potential 
is increased further, the film thickens, which indicates 
a wetting behavior (FIGs. d EH ©• % 

comparing the 

results shown in FIG. [3] and FIG. Q]we readily see that 
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FIG. 4: The local density profiles for water adsorption on 
graphite as a function of reduced distance from the surface, 
z* = z/(T{, at 475 K. The values of x* at which the calcula- 
tions were performed are indicated by the labels in the graph. 
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FIG. 6: The local density profiles for water on graphite as a 
function of reduced distance from the surface at 510 K. Profile 
curves are for X * = 0.731, 0.878, 0.892, 0.919, 0.940, 0.962, 
0.977, 0.995, 1.003, from bottom to top. 
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FIG. 5: The local density profiles for water on graphite as a 
function of reduced distance from the surface at 490 K. The 
values of \* a t which the calculations were performed are 
indicated by the labels in the graph. 



wetting transition of water on graphite occurs somewhere 
between 475 K and 480 K, i. e., T w =475-480 K. At tem- 
peratures in the range of 480 K and 505 K, the prcwetting 
jump in density occurs further to the saturation chemical 
potential with the smaller density jump as temperature 
increases (FIG. [3]). 

At T=510 K, the adsorption isotherm becomes contin- 
uous as the chemical potential increases (FIG.|3|), which 
indicates T pwc <510 K. This is more clearly presented 
by the growth of density profiles shown in FIG. [5] The 
adsorption film builds from a thin to a thick one contin- 
uously with the increase of the chemical potential. The 
density of the first peak in the profiles increases gradually 
to that of an adsorbed liquid. At saturation chemical po- 



tential, the density profile evolves to the one correspond- 
ing to the liquid density except for the first peak adjacent 
to the wall. Comparison of the isotherms obtained at 505 
K and 510 K in FIG. |3] with FIG. 0] indicates that the 
prewetting critical temperature of water on graphite lies 
somewhere between these two values of temperature, i.e. 
T pwc =505-510 K. 

The nature of the prewetting jump of water on graphite 
at temperatures 480- 505 K can be further shown by the 
results obtained from simulations at 490 K, \* =0.992, 
with varying simulation cell dimension in the surface 
normal direction. Shown in FIG. [7] are the density 
profiles obtained from simulations with cell heights of 
H* = h/(if=W, 20, 30, 40, respectively. In order to 
compare the results with consistency, in those four simu- 
lations the area of the graphite wall is kept at 10<7f x lOtTf , 
but the height, or the volume, of the cell varies. It can be 
seen from FIG. [7] that the rapid rise of the film thickness 
to a finite value is independent of the height of the simu- 
lation cell. The film thickness keeps at about 7.5tTf under 
various H* . This is a clear indication that the transition 
is prewetting rather than capillary condensation. We also 
notice that the density profile at H*=10 has more fluctu- 
ations compared with those at H*=20, 30, and 40. That 
could be due to the finite size effects in H* . This confirms 
the importance of using a simulation cell with sufficient 
height in order to obtain reliable wetting transition infor- 
mation. Theoretically the density profiles obtained with 
different H* should coincide with each other. But it is 
hard to achieve this in simulations because of statistical 
fluctuations and metastability nature of the problem. 

One of the clearest demonstrations of the first-order 
nature of many wetting transitions is the observation 
of one stable and one metastable states shown in the 
systems of interest in the temperature range of T w to 
^pwc [IH- The experimental work by Bonn et al. (5fj| 
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FIG. 7: Film density for water on graphite at T=490 K and 
X* =0.992. Curves correspond to varying heights of the sim- 
ulation cell: H*=10, 20, 40, 30, from left to right. 



showed that two different stable values of the film thick- 
ness could be found in the binary liquid mixture of 
methanol/ cyclohexane at a temperature between 295 K 
and 308 K. The thin film of 1.0 nm is the metastablc 
state, and the thick film of 40.0 nm is the stable state. 
Shi et al. [H| observed the similar switching behavior 
between thin and thick films in the simulation study of 
hydrogen isotopes on alkali metal surfaces. 

We also observed the coexistence of stable and 
metastable states for water on graphite, which is shown 
in FIG. [H The probability density function of the wa- 
ter density distribution in the box is collected during the 
simulation with T=490 K, \* =0.968. One major peak 
is located at p* «0.012 corresponding to a thin film, 
at which the simulation samples most frequently. An- 
other small peak exists at p* w0.25, which corresponds 
to a thick film. The bimodal feature of the probability 
density function of the distribution indicates the wetting 
transition at this temperature is first order. The switch- 
ing between the thin and thick films is further confirmed 
by the fluctuation of the number density of adsorbate in 
the box. In the inset of FIG. [8] wc show the evolution 
of water density in the simulation cell as a function of 
simulation steps. At about 11 x 10 6 configurations the 
system abruptly jumps to a higher density of p* «0.25 
from p* «0.012, corresponding to a switching from thin 
film to thick film. This thick film subsequently evapo- 
rates to the thin film at about 14 x 10 6 configurations. 
Comparing the sizes of the two peaks and evolution of 
the number density fluctuation we readily conclude that 
at this state point, the thin film is stable and the thick 
film was metastablc. 

In this work we did not attempt to determine the exact 
values of the wetting temperature or the prewetting crit- 
ical temperature, but only provide estimated ranges for 
them, which are 475-480 K and 505-510 K, respectively. 



One reason is the lack of reliable theoretical method for 
determination of exact T w and T pwc . One of the popular 
techniques being used previously is a power law extrapo- 
lation method [2J] . Theoretical predictions indicate that 
Ap* = (p* vp - p*J ex (T* - T;) 3 / 2 [U. Hence, a plot 
of A/z* versus T* can be used to identify by extrap- 
olating the curve to Ap* = [|, with the saturation 
and wetting transition chemical potentials (/i* vp and p^) 
at various temperatures up to T c of the fluid being cal- 
culated from MHR. However, Shi et al. [3(3 found that 
this power law extrapolation method could be quite in- 
accurate in predicting the wetting temperature of certain 
system. 

Another important concern is the realism of the wa- 
ter potential employed in this study. It is known that 
the transition temperature calculated from simulations 
is sensitively dependent on the solid-fluid interactions 
[H [i^. For example, Shi et al. found that a -10% 
increase in the surface-fluid attraction decreases the wet- 
ting transition temperature of Ar on a CO2 surface by 
3 K |28j . In this work, the choice of the water poten- 
tial will affect the graphite-water interaction implicitly, 
and thus the calculated T w . The SPC/E water potential 
gives by far the best predictions for the critical properties 
of the bulk water among all the available nonpolarizable 
water models, while it still cannot simulate exactly the 
coexistence properties of real water. The accuracy of the 
estimated T w and T pwc may be improved by using more 
accurate polarizable models such as the Gaussian charge 
polarizable model 
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In addition, by modeling the graphite as a smooth sur- 
face, we neglected the possible impact from surface corru- 
gations and dynamics of the surface structure during the 
adsorption. A previous simulation work indicates that 
impact of surface corrugation of the adsorbent on the 
wetting transition behavior of Ne is minimal [25| . But 
it is unclear if the same conclusion is applicable to the 
graphite- water system. 

Cheng et al. [2fJ] proposed a simple model (CCST here- 
after) which interprets the wetting transitions in terms of 
a balance between the surface tension cost of producing a 
thicker film and the energy gain associated with the film's 
interaction with the surface, V(z). The theory results in 
an implicit equation for the wetting temperature 



V(z)dz 



2-) 



Pi - Pv 



(5) 



where pi and p v arc the number densities of the adsorbate 
liquid and vapor at coexistence, 7 is the surface tension of 
the liquid, and z m ; n minimizes the fluid-surface interac- 
tion potential V(z). The wetting transition temperature 
can be calculated by solving the equation since the right 
hand side of Eq. (JS|) is dependent implicitly on tem- 
perature. The impact of solid-fluid interaction on T w is 
reflected by /, and the fluid-fluid interaction is incorpo- 
rated in the model by 7 and pi — p v . 

The wetting transition of water on graphite calculated 
in this work agree quantitatively with the previous pre- 
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FIG. 8: The switching between the thin and thick film of 
water adsorption on graphite at T=490 K, \* =0.968. The 
inset shows the density evolution during the simulation. 



diction [29J using the CCST model, although different 
water potentials are employed in the two studies. It has 
been pointed out by Shi et al. [28| that the wetting be- 
havior predicted theoretically depends on both the well 
depth and well shape of the solid-fluid interacting po- 
tential. Here the potential width is defined as the full 
width at half minimum of the attractive part of the po- 
tential. In Gatica et al.'s work, TIP4P was used instead 
of SPC/E. But we note that the well depth (D) and well 
width (w) of the water/graphite potential for SPC/E and 
TIP4P are almost identical, both with D =9.35 kJ/mol 
and w=0.135 nm, if evaluated at T=475 K and using 
the water dipole moment of 1.85 D. Therefore we expect 
that the wetting transition temperature calculated from 
the CCST model using these two potentials be compara- 
ble. If the simulation results in this work are taken to be 
standard, the CCST model predicted T w of 474 K is very 
accurate indeed. It has been shown that the CCST model 
usually works well in predicting the wetting behavior in- 



volving spherical fluids such as inert gases, but it has 
not been tested extensively with nonsphcrical molecules. 
The fact that this simple model works well in predict- 
ing the wetting of water on graphite, although water has 
a very different kind of potential than inert gases, indi- 
cates that the CCST model contains the essential physics 
of wetting. 



IV. CONCLUSIONS 

In summary, we report the first simulation study of 
the wetting transition of water on graphite surface. The 
wetting transition temperature calculated from GCMC 
simulations is 475-480 K, and the prewetting critical tem- 
perature is 505-510 K. The wetting transition is first or- 
der. The simulation results in this work agrees well with 
the prediction by the CCST model, although the CCST 
model is designed on the basis of the simple fluids such 
as inert gases. 

Finally, we point out that the wetting temperature and 
prewetting critical temperature calculated in this work 
depends on the accuracy of the water potential employed. 
Improvement in the predictions may be made if more 
accurate water potential is available. Future investiga- 
tions can be performed by including the corrugation of 
graphite surface, the finite size effect of the system, and 
by using the more robust simulation techniques such as 
the one proposed by Errington [3l[ . Experimental search 
for the predicted wetting behavior is also warranted. 
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